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We present a hydrodynamical description of the QCD Dirac spectrum at Hnite chemical potential 
as an uncompressible droplet in the complex eigenvalue space. For a large droplet, the fluctuation 
spectrum around the hydrostatic solution is gapped by a longitudinal Coulomb plasmon, and exhibits 
a frictionless odd viscosity. The stochastic relaxation time for the restoration/breaking of chiral 
symmetry is set by twice the plasmon frequency. The leading droplet size correction to the relaxation 
time is fixed by a universal odd viscosity to density ratio rjo/po = (/3 “ 2)/4 for the three Dyson 
ensembles /3 = 1,2,4. 

PACS numbers: 12.38Aw, 12.38Mh, Tl.lOPm 


1. Introduction. QCD breaks spontaneously chiral 
symmetry with the emergence of an octet of light mesons 
that permeate most of the hadronic processes at low en¬ 
ergies [ij. Dedicated lattice simulations are now in full 
support of this spontaneous breaking [^. Fundamental 
light quarks become constitutive and heavy producing 
most of the mass of the elements around us. 

A remarkable feature of the spontaneous breaking of 
chiral symmetry is the large accumulation of the eigen¬ 
values of the Dirac operator near zero-virtuality with the 
formation of a finite vacuum chiral condensate Q . Small 
eigenvalue virtuality translates to large proper time, as 
light quarks travel very long in proper time and delocal¬ 
ize. The zero virtualiW regime is ergodic, and its neigh¬ 
borhood is diffusive [j. This behavior is analogous to 
disordered electrons in mesoscopic systems 

The ergodic regime of the QCD Dirac spectrum with 
its universal spectral oscillation is described by a chi¬ 
ral random matrix model In short, the model sim¬ 
plifies the Dirac spectrum to its zero-mode-zone (ZMZ). 
The Dirac matrix is composed of hopping between N-zero 
modes and N-anti-zero modes because of chirality, which 
are sampled from gaussian ensembles thanks to the cen¬ 
tral limit theorem. The model was initially suggested as 
a null dynamical limit of the instanton liquid model 0. 

QCD at finite chemical potential /r is notoriously diffi¬ 
cult to sample on a lattice due to the sign problem Q . A 
number of chiral models have been proposed to describe 
the effects of matter in QCD with light quarks [l|. In vac¬ 
uum, the chiral random matrix model simplifies the QCD 
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Dirac spectrum to its ZMZ. In matter, the light quark 
zero modes are involved. Their chiral and cross-hopping 
in the ZMZ is suppressed exponentially, and the corre¬ 
sponding Dirac matrix is banded and not random. How¬ 
ever, large matter effects reduce the banded matrix to its 
diagonal, localizing the quark zero modes into molecules. 
In the I-matrix model the chiral random ensemble is de¬ 
formed W a constant matrix, leading to localization at 
large p [l3 ■ In. the 2-matrix model the deformation 
is still random and only generic for moderate /r with no 
strict banding at large p, [111 . Il2| . The 1-matrix a ppr oach 
to QCD at finite fi has been discussed by many [ll. flilLIl4| . 

In this letter we would like to combine the ergodic char¬ 
acter of the chiral random matrix model for the low-lying 
modes in the ZMZ with the universal character of the hy¬ 
drodynamics approach, to describe the relaxation of the 
QCD Dirac eigenvalues at finite p.. We will obtain the fol¬ 
lowing new results: 1/ a hydrodynamical description of 
the Dirac eigenvalues as a droplet in the complex 2-plane; 
2/ small amplitude deformations in the droplet that are 
gapped by the emergence of a plasmon with an odd vis¬ 
cosity; 3/ a non-perturbative estimate of the stochastic 
relaxation time for breaking/restoring chiral symmetry 
in matter. 

2. The model. A useful model for the eigenvalues 
of the QCD Dirac operator at finite /r makes use of a 
2-matrix model [ll|, [l^ 

. N N 

Zii[mf]= J - z'jf 

i=l i<j 

x{zj + m))^i ( 1 ) 

for quarks in the complex representation or /? = 2 with 
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a = /3(iy + 1) — 1. ly accounts for the difference between 
the number of zero modes and anti-zero modes. The 
potential is 


W(z) = ^(lzl^-^(z^ + z^) ( 2 ) 

with P = 1 — T = 2/r^/(l -I- fi^). For ^ 0, r Ri 1 

and P Ri 2/r^, so that fF(z) ~ —{N/^‘^){z — z)^, which 
restricts the eigenvalues to the real axis. 

For the other quark representations with /3 = 1,4 the 
joint distribution is more subtle [l^. Throughout, (ED 
will be assumed for /? = 2, but all results extend to /3 = 
1,2,4 for large N and/or the quenched approximation. 
All units are expressed with a = 1 unless noted, which 
is related to the vacuum chiral condensate through the 
Banks-Casher formula 0. Specifically ^/a = Ig'^qlo/n = 
1, with n = N/V 4 the density of zero modes. 

ED can be re-written as an average of the complex 
fermion determinant 


N 

^/3[^7^/] = / Y[d^Zi{zf l’®'o[2]P (3) 

i —1 

using the real many-body wave-function which 

is the zero-mode solution to the Shrodinger equation 
= 0 with the self-adjoint Hamiltonian 


1 ^ 
i —1 


■ a,. 


( 4 ) 


Here di = d/dzi and the gauge potential is ai = diS with 
S'[z] = — In^'njz]. In (j?]) the mass parameter is to = 1/2. 

Following [l6|, we observe that the Vandermond de¬ 
terminant A = ni<j induces a diverging 2-body 

part in Hq. Using a similarity transformation, we can 
re-absorb it in 4/ = 'ko/'\/A, and the new many-body 
Hamiltonian is 


H=^HoVA (5) 

We will refer to Q as the quenched Hamiltonian. The 
phase quenched Hamiltonian follows a similar reasoning 
by rewriting m as 


r N /z2 , 2 \ 2 

Zp[mf]= J Y[d^z, - 2^^2 ) ( 6 ) 

Below we detail the hydrodynamical construction in the 
quenched approximation, and quote the minimal changes 
in the phase quenched approximation for /3 = 2. 


re-write m in terms of the density of eigenvalues as a 
collective variable p{z) = ~ ^i)- After some 

algebra and up to boundary and ultra-local terms, we 
obtain 

H = Jd^z p{z) ^ ((V7r)2 + (A)2) = J d^zh (7) 

with the pair tt, p canonically conjugate. We will restrict 
our discussion to the classical limit with the pair obeying 
the Poisson brackets {tt{z), p{z')} = (5^(z — z'). Defining 
the even density p^{z) = p{z) + p{—z), we have 

A = A+iv(/3p^(z) + (/3-2)lnv^) (8) 

Here pL is the logarithmic transform of p 

[p]l = Pl{z) = J dz'\n\z - z'\p{z') (9) 

and the vector potential (t± = 1 ± t) 

7* NB, . a , , 

A=--^(r-x,r+y) + 2|^(a:,?/) (10) 

The equation of motion for p yields the current con¬ 
servation law and the Euler equation for v. Defining 
mv = Vtt, they are specifically given by 


dtp + V ■ {pv) = 0 

dtTT + -mv + — 

2 2m 

(pA) --^V-[p^Ali = 0 (11) 

The steady state flow from EID corresponds to Bernoulli 
law with dtTT = C a fixed constant. Note that all the 
relations hold for large but finite N, provided that the 
fluid density and velocity are sufficiently smooth. 

4. Hydrostatic. The quenched hydrostatic solution 
is encoded in the condition A(z) = 0 and tt = 0. Using 
the formal identity pL = (27r/V^) p, we have 


where the integration constant k = 1 + a/{ 2 NjB) is fixed 
by the density in leading order, and A is the area of the 
eigenvalue density. 

In the phase quenched approximation for /3 = 2, the 
vector potential EHD is now shifted 


3. Hydrodynamic. In the quenched approximation, 
we can use the collective coordinate method in to 


^ Mr ^ 

A A+ -;^Vln|z^ + 


m}\ 


(13) 
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with the hydrodynamical equations m unchanged. The 
corresponding phase quenched hydrostatic density m is 
modified in subleading order 


p{z)^p{z)-^VHn\z^+m}\ (14) 

Both the quenched and phase quenched approximations 
describe an elliptic droplet at large N. 

5. Droplet boundary. It is useful to recast the 
hydrostatic equations in complex form at large N 


z 


z 




(15) 


S{z) is a Schwartz function with S{S{z)) = z dl, i.e. 
S{z) is the inverse of S{z). S{z) is analytic outside A, 
with S{z) fv TZ + P/z asymptotically. We note the simi¬ 
larity with the Blue's function and its use in determining 
the domain boundary through conformal mapping in M- 
The support of the eigenvalues A is determined by a 
conic curve 


\z\^+ C{z'^+i^)+c'{z + z)+c" = Q (16) 

Using the substitution z = S{z) with the asymptotic 
form of S{z) fixes the constants in (I16|) . The domain in 
z = (a;, y) is an ellipse ja^ + y'^IaP_ = 1, with the axes 


a\. 1 ± r 
2P 1 =F T 


(17) 


The ellipse remains un-split with A = 7ra+a_ = 2'kP. 
The area is preserved under r-deformation in ([2]). 


6. Plasmons. It is useful to analyze the small defor¬ 
mations in the density and velocity profile by linearizing 
the current conservation law in (in]),i-e. dtSp+poS/'^STT = 
0, which is readily solved using 


with 


W[4>] 


V 


/3-2 (5p\ 

4 Po/ 


( 21 ) 


Using again the formal identity /l = (27r/V^)/ and 
defining the small longitudinal field (^= we obtain 


N / d?zdt-^ 
J 2m 

X 


( 22 ) 


t^PPo -.v , /3 - 2_2 - 


after the rescaling Nt —>■ t. The small longitudinal 
excitations in p are gapped by the plasmon frequency 
ojp = 2Tr(3po/N. For an elliptic droplet of large area A, 
(1^ leads to the quadratic dispersion law 


uj{k) Ri ± 


UJp 



(23) 


Here |fc| is conjugate to |z|. The gapped spectrum means 
that the droplet is uncompressible. For /3 = 1,2 with 
quarks in the real and complex representation the branch 
(l23l) describes a plasma fluid. For /3 = 4 with quarks in 
the quaternion representation, (PI shows the start of a 
roton-like branch a possible indication of superfluidity. 
For that the higher order k'^ term is needed. 


7. Odd viscosity. There is an interesting analogy 
between the droplet of Dirac eigenvalues at finite chem¬ 
ical potential, and the quantum Hall effect as a fluid of 
neutralized charged electrons in the plane [l^ To 
illustrate the analogy, we first note that sources the 
magnetic field 


B{z) = V X A* = — - 7ra(5^(z) (24) 

with the dual notation V* = f-ijVj subsumed. Amusingly, 
(HI describes a Coulomb fluid in the magnetic field (1^ . 
In large N the density of eigenvalues is uniform 


5p = -poV‘^4> 5 tt = dtc, 
Inserting (fTSll in the canonical action 

S = y cfzdt{Trdtp—h) 
yields in the quadratic approximation 




(18) 


Piz) 


N 

2ttP 


izB 

27r 


(25) 


(19) 


( 20 ) 


which is the density of a quantum Hall droplet with filling 
fraction v = 1//3. The plasmon frequency is the cyclotron 
frequency ujp = B/M with M = N the analogue of the 
effective mass. I identifies with the magnetic length. 

The ^^-contribution in (1^51) is reminiscent of the odd 
viscosity in the fractional quantum Hall effect. To show 
this, let fr = in and define the collective velocity mv = 
Vtt-I-A*, then (I7| is a free flow-like Hamiltonian modulo 
ultra-local terms. 
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J (f z p{z) ^ v~^ ■ V (26) 

In our case ^ v but in the fractional quantum Hall 
effect they are equal, making A* a real gauge-field and v 
a real and gau ge-i nvariant flow velocity for flux-riding 
quasi-particles |l9l . [ 23 . Assuming that n and p are 
canonical and after some algebra, the Euler equation fol¬ 
lowing from (1261) yields the momentum conservation law 


dtipmvi) + VjTij = 0 (27) 


with the stress tensor 


We note that similar non-linear equations emerge from 
the diffusion of non-hermitean matrices [ 2 ^. 

The solution to (1321) with a free boundary or large 
droplet size A can be obtained using the method of char¬ 
acteristics. Specifically 


dts 

ds 

dr 

ds 

ds 


-r{l3f 

fW- 


Nj5r 
N fir 


(33) 


with the conditions t£:(s = 0) = 0, r{s = 0) = rg and 
/(s = 0) = /(ro). For f{r,tE =0) = Kr with K » po, 


8 — 2 

Tij = mpViVj -\ - - — p (ViW* -I- V*Vj) (28) 

The first contribution is the classical free fluid part. The 
second contribution is the odd viscosity contribution fol¬ 
lowing from the breaking of parity in 2-dimensions [2lj, 
with 


Vo_ ^ /3- 2 
p 4 



(29) 


tE = —as 


N8 


In 


ro -b a - (ro - a)e"i^ 


2a 


ro -b a -b (ro - a)e i 
r = a-« 

ro -b a - (ro - a)e i 
a=^2Krll^/N 




(34) 


which is the coefficient of the term in (l23l) . In the 
quantum Hall fluid, po originates from a mixed gauge- 
gravitational anomaly We note that the pair v,v* 
are orthogonal. This explains that the fc^-contribution in 
(l23l) acts as the even (shear) viscosity but without the i 
for dissipation. No vorticity is therefore expected. 

8. Instanton and relaxation time. We identify 
the zero energy configuration in ([T]) as an instanton so¬ 
lution with imaginary velocity with tt —>■ in, i.e. 

iVTTp - |A|2 = 0 (30) 

that satisfies the analytically continued in time conser¬ 
vation law (t —>■ —He) 


We first use the three equations in (IM)) and solve for ro = 
ro(r, tfi). We then substitute the answer in the fourth 
equation in (IM)) to find explicitly /(r, t^;) in general. Its 
large time asymptotic for s —>■ —c» is Ie ~ —as and 



(1551) relaxes as to /(r, 00 ) = Nr/2P lead¬ 

ing to the hydrostatic density po = N/{2nP). We iden¬ 
tify Tr PS l/2u!p with the relaxation time after rescaling 
NIe —^ tE- The scalar density relaxes by emitting two 
longitudinal plasmons. In physical units 


-dt^p + V ■ {pVn)=0 (31) 

Without loss of generality and for simplicity we choose 
T = 0 in ([5]) so that the hydrostatic droplet is circular. To 
solve m we set p(0, z) = K/n ^ Po, which corresponds 
to all eigenvalues localized in a small disc centered around 
the origin to simulate a chirally restored phase at finite 
p. dSH) simplifies by radial symmetry 


drPL{r,tE) = f{r,tE) 

rdtj^f + r{l3f - ^^)drf + fiPf 


(32) 


N8r 


) = 0 


after re-instating 1 = y/a = |a’^ a|n/ n, and adding the 1 
to reproduce the p = 0 result in [1^ . A simple extension 
to finite temperature amounts to a re-definition of units 
or y/a —^ y/Er = \q^q\T/nT as in [iElil. 

Finally, we note that Tr Ri l/2a;p translates to a dif¬ 
fusive time with P = Al2n ps 2/3Tr. The diffusion 
constant is D = 2/3. An estimate of the finite droplet 
size corrections follow from (1551) using the substitution 
ojp ^ uj{k PS 1/-\/A). The leading correction is controlled 
by the odd viscosity to density ratio in ((55)) and is small. 


2/2 
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9. Conclusions. The hydrodynamical description 
captures some key aspects of the QCD Dirac eigenval¬ 
ues in the diffusive regime at finite chemical potential. 
It supports an instanton that describes the stochastic re¬ 
laxation of the Dirac eigenvalues as a fluid. The fluid 
is uncompressible and exhibits dispersive plasmon waves 
that can be used to estimate the time it takes for a chi- 
rally symmetric phase to relax to a chirally broken phase 
in matter. The time estimate is non-perturbative and 


gauge-independent. 
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